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Abstract 



In this paper we attempt to define Grassmann numbers as elements living in multi- 
dimensional space, that exist independently of path integral. Furthermore, we define 
path integral as literally a "limit of the sum" in such a way that it approximates the 
expected mathematical properties of Grassmann integration. We then proceed to show 
that "Harmonic oscillator" in Fermionic space produces a linear function. The latter 
has two degrees of freedom per variable (the slope and the height at the origin), just 
as expected from Fermi exclusion principle. At the same time, the linear function is 
living in a continuum space; in other words, only two out of uncountably many degrees 
of freedom are "utilized". However, we argue that we would be able to utilize the 
rest of the degrees of freedom if we were to modify the structure of our Lagrangian. 
Such modification involves the use of functions that can not be generated from the 
usual wedge products of Grassmann numbers. Thus, it requires an appleal to the 
fact that Grassmann numbers are well defined mathematical objects, independently of 
standardly-accepted formalism. 



1. Introduction 



It is well known that, in Feynmann path integral picture, Grassmann numbers are used in 
order to define fermionic field. At the same time, the integral of Grassann numbers is not 
viewed as a limit of a sum. Instead, it is believed to be a "formal" operation of symbols. 
This is probably the case for two reasons. From the mathematical point of view, as will be 
outlined in the next section, naive attempts to view Grassmann integral as literally a "limit 
of a sum" would create some paradoxes that would need to be addressed. On the physical 
side, Fermi exclusion principle predicts only two states per variable. Thus, if we were to 
assign to it a continuum spectrum (as would be required in order to define integral in a 
"usual" way) we would have uncountably many degrees of freedom, out of which only two 
will be utilized. At the same time, however, it is common place to "pretend" that Grassmann 
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variables are somehow "physical". Examples of this range from the "rotations" of fermionic 
fields in standard model, to the physics beyond standard model, including supersymmetry. 

Of course, there are ways to justify the "conventional" outlook on Grassmann numbers. 
For example, one can point out that in Fock space formulation the creation and annihila- 
tion operators are equally well defined for bosons and fermions, and neither require path 
integral. Furthermore, it is possible to to redo the entire supersymmetry without reference 
to Grassmann numbers altogether ([2] and [1]). Now, it is logically possible to continue to 
use Grassmann numbers and path integrals we "dont believe to be real" and at the same 
time continue to hold true to realistic interpretation of quantum mechanics. After all, when 
we perform an addition by writing numbers in a column and drawing a line in the bottom, 
we don't provide physical meaning for either the position on the column or that line. This 
is merely a tool to compute the sum; yet, we do know physical meaning of the latter. In 
principle it is possible for a realist to understand, and use, Grassmann numbers and their 
"integrals" in that same way. The Fock space formulation is "analogous" to the "sum", 
which we do know the meaning of, while path integral is " analogous" to writing numbers in 
column which has no meaning. 

Nevertheless, while the above allows us to continue to be realist, it imposes some lim- 
itations on the possible "realistic" viewpoints we can consider. As we have just stated, we 
were forced to view Fock space as more fundamental than path integral. In curved space- 
time, however, this would imply "preferred" foliation of spacetime into hypersurfaces and 
"preferred" definition of particle number, contrary to the spirit of Bogolubov's transforma- 
tions. We can then attempt to justify the preferred frame by claiming that we "happened" 
to believe in one anyway for the "independent" purposes related to quantum measurement 
problem. The other alternative would be to claim that the reality is to be computted through 
Feynmann diagrams whose in- and out- states are bosonic, while fermions only come into 
the loop terms. We can then re-interpret Grassmannian integral along "fermionic line" as a 
"two-point function" on bosonic line, the Grassmannian integral on fermionic loop as "one 
point function" on bosonic line, etc. which would be defined independently of fermionic 
fields. In the context of pilot wave models, this thinking was used by Struyve and West- 
mann [9]). While it is possible to argue that fermions don't exist and we "falsely infer" 
their existance by looking at bosons they emit, this imposes some limitations on the theories 
we can consider. For example, Valentini in his thesis attempted to use fermionic fields as 
beables which would now have to be deemed invalid. This begs the question: do we really 
want to continue to impose restrictions similar to above, or is it perhaps possible to realisti- 
cally define Grassmann numbers which would make us free to consider broader spectrum of 
theories? 

In fact, the argument that "everything involving Grassmann numbers can be explained 
away by other means" can be turned around. After all, we don't have "direct evidence" for 
anything else either. Instead, we gather a pattern of "coincidents" and, once that pattern 
is long enough, we conjecture that whatever this set of "coincidences" point to, indeed, 
exist. From this point of view, the "coincidence" that Grassmann numbers provide excellent 
"computational tools" for topics as diverse as path integral formulation of standard model 
as well as supersymmetry seem to suggest that perhaps they do, in fact, exist as literal 
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mathematical object. One thing that draws a big hole in this evidence is the fact that we 
would not expect the behavior of "limit of the sum" to match the Grassmann integral as 
we know it. This is the main obstacle we overcame in this paper. Indeed, we have shown 
in Section 3 that if we assume two different kinds of products we will specify, and if we are 
to perform an integral over a contour that resembles "random walk" in a space with very 
large number of anticommuting dimensions, and, finally, if we assume a metric as defined 
in Equation [13 we will, in fact, predict that "limit of the sum" will obey the Grassmann 
integration properties! 

As mentioned earlier, the idea of defining Grassmann numbers would probably meet 
a physics objection that their "continuous" spectrum would contain a lot more degrees of 
freedom than could be directly observed. At the same time, however, it is possible to claim 
that from pure kinematical point of view we could observe these degrees of freedom, after 
all. We will claim, however, that the fermionic "dynamics" predict that the system with 
infinitely many degrees of freedom will immediately "collapse" into the linear combination 
of only two of them, and stay there. This will be shown to be a consequnce of Grassmannian 
"path integral" taken over the "continuous" anticommuting space we have defined. At the 
same time, this will, in principle, allow for the possibility that under some modifications of 
the theory (in particular, the way in which we treat "time derivative" term) such collapse 
would either be slowed down or temporarily reversed, thus leading to non-trivial physical 
predictions. This idea might particularly be fruitful in case of quantum measurement. 

2. Key mathematical difficulties and how to avoid them 

Before we proceed to define Grassmann numbers, let us address some of the mathematical 
arguments as to why the pursuit might seem "impossible". In other words, let us first 
convince the reader that what we will do does not imply any mathematical contradiction, 
and after that, in the subsequent sections, we will actually do it. Let us, therefore, make a 
list of apparent "contradictions" we would face (a through f) and after that we will "answer" 
each of the "contradictions" in the corresponding items a' through f: 

a) Anticommutativity implies that the product of any two co-linear elements is zero. 
Now, if a one- variable integral is over a line, then multiple integral is over all possible choices 
of n points on that line. Since these points are bound to be colinear, their wedge product 
will return zero, contrary to what we want. 

b) Even though products of two Grassmann numbers commute, (9i9 2 )9 3 = +9 3 (9i9 2 ), 
they can not be identified with real line. After all, {9\9 2 ) 2 = 0. Thus, it is a bit puzzling 
why the integration of "products" 9d9 would, in fact, return a "real number" (namely 1). 

c) In the definition of the integral we demand that the integral over odd function is 1 
while integral over the constant is 0. On the other hand, if we were to view integral as a 
sum, we would be expecting just the opposite. 
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d) The fact that 

J d9 1 d9 2 9 1 9 2 = +1 (1) 

with plus rather than minus on the right hand side, implies that we assumed that (d9id9 2 )(9i9 2 ) = 
+ (d9i9i)(d9 2 9 2 ), contrary to what anticommutativity would imply. 

e) In order to define the integral as "limit of the sum", we need the notion of "norm" 
since the latter is part of the definition of the limit. But, if we define norm in a "naive" way, 
\9\ 2 = 99, then anticommutativity will imply that \9\ = for all 9. This would render the 
definition of limit and, therefore, integral, meaningless. 

f) Scaling symmetry implies that the integral of e a9 is a _1 e a0 ; yet, Grassmannian integral 
is ae ad instead. 

Because of these difficulties, it is commonly assumed that Grassmann numbers are not 
mathematical objects at all; they are merely "symbols", and their integral is likewise viewed 
to be "manipulation of symbols" as opposed to "limit of the sum". In this paper, however, 
we found a way to get around these difficulties and ultimately come up with "realistic" view 
of both Grassmann numbers and their integrals. In particular, we adress the above list of 
issues in the following way: 

a') Instead of assuming that one- variable integral is over a line, we assume that it is over 
a contour that is living in infinitely dimensional space, spanned by {ei, e 2 , ■ • • }■ The contour 
has some specific afore-given, non-trivial shape, which is expressly different from the straight 
line. The n-variable integral will be the integration over n random points on the contour. 
In light of infinite dimensionality, combined with the curvature of the contour, the vectors 
pointing towards these points are linear independent. Thus, their product is no longer zero. 

b') We have two kinds of products: dot product (•) and wedge product (A), where we 
assume that only the latter anticommutes, but not the former. The typical expression under 
the integral takes the form d9 1 ■ d9 2 ■ (9i A 9 2 ). The presence of dot-product prevents us from 
saying that wedging it with itself will return zero. Thus it is conceivable that the above 
product is, indeed, a real number. 

c') We can turn the integral over even function into and over odd function into 1 by 
using one of the two tricks: 

(i) We can displace the center of the region of integration away from the origin. As a 
result, both even and odd functions will integrate to non-zero values. After that we can 
introduce a sequence of these regions and take the limit over the corresponding sequence of 
integrals. It is conceivable to select the sequence in such a way that the limit will return 
zero for even functions but not the odd ones. 

(ii) We can, contrary to part (i), assume that the center is still at the origin. However, 
we can impose non-trivial measure that changes sign; in fact, the measure that we impose 
happens to be an odd function! As a result, the product of measure with the odd and 
even function will be even and odd respectively. This would imply the desired result: even 
functions would integrate to zero while odd ones would not. 
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d') This is merely the issue of convention. We can go back between the two conventions 
by defining an extra product operation, o, by Xi ° Xi — Xi " Xi- Thus, in light of transitivity 
of •, we have d9 1 ■ d6 2 • (6*i A 9 2 ) — (d9 2 o dOx) ■ (6>i A 6 2 ). We will agree that 

1 = J (d0 1 o dd 2 ) ■ (6 1 A 9 2 ) = J d0 2 ■ d0! ■ (9 1 A 2 ) (2) 



e') We will introduce a definition of norm independant of products. 

f) As stated in a', we are integrating over a contour that has non-trivial shape. The 
shape of the contour violates scaling symmetry. Now, in order to prove that integral of e ad 
is a~ 1 e ad , one needs to appeal to scaling symmetry. Since the scaling symmetry is no longer 
valid, neither is alleged property of the integral. 

In light of the fact that we have gotten rid of "no go" arguments, we will now attempt 
to go ahead and define Grassmann numbers as literal mathematical objects that exist inde- 
pendently of integration, and which might have a wide spectrum of continuous values. We 
will also define integral as a "limit of the sum" and demonstrate that, if the dot product 
is defined in a certain way, the outcomes of the integrals over contours that make conse- 
quetive terms in different dimensions (inside almost- infintie dimensional space) can in fact 
be expected to return the typical properties of Grassmann integral. We will then work out 
the harmonic oscillator in this formulation. Finally, in Conclusion we will briefly mention 
possible physics applications of our "new" definition of Grassmannian integral, that are to 
be explored in future research. 



3. Grassmann Integral 

As was mentioned in part a', anticommutativity is accomodated through a very large dimen- 
sionality. In particular, our space will be spanned by unit vectors e±, ■ • ■ , ejv. Furthermore, 
as was mentioned in part b', we are going to introduce two kinds of products: dot-product 
and wedge-product. The two products coincide when it comes to orthogonal unit vectors: 

i 3 =>■ e i • e j = e i A e j (3) 

On the other hand, if we are to multiply the unit vector by itself then the wedge product 
will give us zero, whereas the dot product will give us 1: 

e k A e k = , e k ■ e k = 1 (4) 

It is easy to see that the same object has more than one dot-product representation (for 
example, 1 = e 5 • e 5 = e 7 • e 7 ). At the same time, wedge product representation is unique 
as long as we make sure to write it down in an order of increasing of indexes (thus avoiding 
the confusion between e\ A e 2 and e 2 A ei). From the above properties one can compute dot 
products of more complicated objects. For example, 

(ei A e 5 A e 7 ) ■ (e 4 A e 5 ) = {-e 1 A e 7 A e 5 ) • (-e 5 A e 4 ) = {e 1 A e 7 A e 5 ) • (e 5 A e 4 ) = (5) 



5 



= (e 1 -e 7 -e 5 )-(e 5 -e4) = (ei-e 7 )-(e 5 -e 5 )-e4 = (ei-e 7 )-l-e 4 = ei-e 7 -e 4 = eiAe 7 Ae 4 = — eiAe 4 Ae 7 

It should be noted that the consistency of the procedure such as above needs to be formally 
proven. Strictly speaking, we have to define dot- and wedge- products of the "chains" of 
arbitrary length and then prove their associativity. Such proof, however, turns out to be 
surprisingly long and for that reason we have decided to leave it out. A Grassmann number 
9 is given as a linear combination of unit vectors without any products: 

= Y^ X ke k (6) 

where \ k are commutting reals while e k anticommute as just described. When we are per- 
forming Grassmann integral, we are integrating over afore-given contour C as opposed to 
entire domain. That contour is given by some agreed upon set of differential functions 
Afc(cr): 

9(a) = J2^(o-)e k (7) 

Now we will assume that we have very large, but finite, dimensionality, N. At any given 
time, only n out of N coordinates evolve, where n is some small integer (such as 2 or 5). For 
example, if n — 4 then our contour will start out parallel to X1X2X3 hyperplane, then it would 
slowly "rotate" to the direction parallel to £2X3X4 hyperplane, then it would move parallel to 
X3X4X5 hyperplane, and so forth (this agrees with the "non-trivial shape" we mentioned in 
part a'). During the time of "shift" from hyperplane to all four coordinates 

evolve, which is why in this particular example we say n = 4 rather than n — 3. We will, 
therefore, assume that the total increment of each coordinate is reasonably small 

\ k (a = 00) — Afc(<7 = —00) = AA (8) 

Now, if we are perform the integral of d9 over the above contour, we obtain 

N 

/^ = AAVe fc (9) 

On the other hand, the integral of d9 ■ 9 is given by 

J^d9-e = ^2(e i -e j J ^Xj(a)d<Tj (10) 

Now, the above sum goes over both i ^ j and i = j. The i ^ j terms will produce 
ej • ej = et A ej, while i = j terms will produce e k • e k — 1. Thus, we obtain 

Jd9.9 = J2(f ^(^) + E ( e * A ^ J ^ A >)^) ( n ) 

By evaluating the first integral on the right hand side, we obtain 

fiS.B = + yj („ A e, f d ±\W)da) (12) 
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Now, if we take " Pythagorean" absolute value of Equation [91 we obtain 



Pythagorean 



de 



iV(AA) 2 (13) 



which happens to be twice larger than the integer part of Equation [121 In light of the fact 
that we "want" the left hand side of Equation [9] to be zero and the left hand side of Equation 
12] to be 1, this doesn't appear too optimistic. However, situation improves drastically if we 
replace the Pythagorean metric with "Special Manhattan" one (where "special Manhatten" 
differs from Manhatten in that we take maximum of coordinate values instead of summing 
them): 



Pythagorean =^ | ^ ^ e t| = y X l ( 14 ) 
Manhatten =$>■ | \k e k\ = A& (15) 

Special Manhatten =>- | ^^Afcefej = max{Afc} (16) 

If we now adopt the " Special Manhatten" metric in our Grassmann space, and generalize it 
to tensors according to 

| Yl A «i-«fe e «i A ■ ■ " A e ^ = max{A 01) ... )( j£; e N} (17) 

k ai,---,afe 

we will, indeed, conclude that the real part of right hand side of Equation[9]is "much smaller" 
than the real part of right hand side of Equation [T21 as desired. The only thing that is left 
to do is to select A A in such a way that the latter is equal to 1. One can easilly see that this 
is accomplished through 

AA=y| (18) 

Let us now look at the integral of two Grassmann variables. We would like to introduce 
two versions of double integral, 

I d6 l ■ dd 2 ■ (0! A 6 2 ) = ( e a ■ £b ■ (e c A e d ) I d ( T 1 d ( T 2 A , a ( ( T 1 )A^(a 2 )A c ( ( 7 1 )A d (a 2 )) (19) 



abed 



and 



/ d9 l ■ d9 2 ■ {9 l A9 2 ) = J2 (e« • e b ■ (e c A e d ) f rf ( x lC / 2 A , a ( ( T 1 )A^(a 2 )A c ((T 1 )A d (a 2 )) (20) 



abed 

Now we would like to represent both integrals as 



de 1 • 0ij (J de 2 ■ e 2 j = (21) 

e b ■ e d / da 2 X' b (a 2 )X d (a 2 ) 
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This would require us to assume 



3fJ(e a • e b • (e c A e d )) = WANT K(e a • e c )K(e 6 • e d ) (22) 

(e a A e 6 ) • (e c A e d ) = WANT K(e a • e c )K(e 6 • e d ) (23) 

where 3? stands for "real part". The above, however, is not correct. For exaple, if we set 
a = l, b = 2, c = 2 and d — 1, we obtain 

»(ei • e 2 • (ea A ei)) = »((ei A e 2 ) • (e 2 A e x )) = »(-l) = -1 (24) 

and, at the same time, 

3?(ei • e 2 )9R(e 2 • e x ) = x = (25) 

Intuitively, what we are facing is "unwanted contraction". We would like to separately 
"contract" dQ\ with 9i and d6 2 with 9 2 , leading to separation of variables. But, unfortunately, 
we obtain "unwanted" contractions of dQ\ with 9 2 and d6 2 with d\ in the process. These 
contractions are also real-valued and, therefore, special Manhatten metric is no longer of 
any use. Instead, we get unwanted large real component in our integral. 

However, we still have one "trick" left! In particular, we can design the shape of the 
contour in such a clever way that the "unwanted contractions" occur only when <j\ and 
a 2 are relatively close to each other. As a result, the part of the integral "without" such 
contractions will simply "outweigh" the part of the integral "with" them. In particular, we 
will postulate 

dXk , „ /, n\ /, rr 



da*°^{ k -2h <Cr< l k+ 2h (26) 
Now, if we assume that the contour starts at the origin, 

A(-oo) = (27) 

this will imply that for i <C j we will have d6i ■ 8j contraction, but we will not have ddj ■ 9{ 
one. Now, in the exaple of Equations [2H and [25j we need both contractions present at the 
same time in order to produce a real number. On the other hand, if one is present while 
the other is not, it is easy to see that we will get some "flowing" components, whose total 
contribution will end up being "very small" thanks to special Manhatten metric. 

What we have just stated about two-variable integral can be generalized to several 
variable case. However, that generalization relies on the fact that — o~j\ ^> na throughout 
the "dominant" part of o\ X • • • x o~d (after all, this is a prerequisite for the claim that, per 
Equation [261 there are no unwanted contractions). The latter, in turn, holds true only if 

D < iV (28) 

where D is the number of variables we are integrating over, 

J de 1 -...-de D Af(e 1 ,--- ,e D ) (29) 
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while N is the number of variables the contour lives in: 

6 = X 1 e 1 + --- + X N e N (30) 

Geometrically, we are integrating over arbitrary distribution of D points on a contour that 
makes the order of iV "turns". If two of these points we integrate over accidentally falls 
within the same "turn" then we would have unwanted contraction. If they fall within two 
neighboring turns, we might still have unwanted contraction due to the "directions" of these 
"turns" sharing some of the coordinates but not all; in fact, this will be a concern even if 
they are distant neighbors, as long as they are either n-th neighbors of each other or closer. 
If, on the other hand, the points are placed at " far away" turns, the tangent lines towards 
the curve at these points do not share any of the coordinates and, therefore, there are no 
unwanted contractions. If D <C iV then the latter scenario is almost guaranteed to happen. 
If, on the other hand, N <^ D, then we are forced to face the former scenario and deal with 
unwanted contractions. 

Now, in the real physical situation we have separate fermionic coordinate associated with 
every point in space. Thus, in order to say that D < iV we have to first assume that the 
spacetime is discretized and then further assume that the number of discrete points is much 
smaller than N. These two assumptions are quite separate. As far as discretization itself, it 
has to be assumed in bosonic case as well since integral with infinitely many dimensions is 
simply not defined. But the second assumption that the number of lattice points is much 
smaller than iV is strictly a consequence of Grassmannian nature of the integral. Of course, 
since we have not specified N, any "finite" number of lattice points will have to be "much 
smaller" than something. Thus, the statement that we have made might be interpretted as 
simply saying that " N is bounded below" as opposed to " D is bounded above" . The fermion 
doubling issue, of course, needs to be addressed; this, however, is beyond the scope of this 
paper. 

4. Complex extension and linear combinations 

As the reader probably noticed, we have defined Grassmann integral as the integral over 
the contour rather than the entire domain. This raises an important question: how can 
we "rotate" Grassmann numbers? Yet, rotation of Grassmann numbers seems to be a very 
important property. First of all, it is common place to perform integral over 6 = 9± + iQi 
and 8 = 0\ — i9 2 , instead of directly integrating over 9\ and 82- Even more importantly, in 
gauge theory, the symmetry transformations are commonly applied to fermions; and they 
would not be defined if we couldn't rotate Grassmann numbers. We propose to address the 
above by making clear distinction between Grassmann numbers and Grassmann variables as 
follows: 

a) Grassmann numbers are defined in the entire N- dimensional space, and are not con- 
strained to any contour nor any other subspace of that space. Likewise, all of our products 
and sums are defined on the space of all Grassmann numbers. 

b) The term variable is used only if a Grassmann number under discussion is constrained 
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to some agreed upon contour. 



c) We can, in principle constrain different Grassmann variables to different contours 
from each other; or we can choose to constrain them to the same contour. It is up to us! 

d) A given "variable" (constrained to afore-specified contour)' is said to be Grassmann 
variable if the contour it is constrained to leads to the integration properties described in 
the previous sections, up to the agreed-upon approximation. 

e) A variable that violates part "d" by living in the "wrong contour" can still be referred 
to as "variable", as long as we are not calling it "Grassmann variable". Likewise, we can 
still integrate over that "wrong" contour, as long as we don't expect the results to match 
with Grassmann integration. 

f ) Since Grassmann numbers are defined independently of the presence of any contours, 
the "variables" that violate part "d" by living on the "wrong contour" are still Grassmann 
numbers even though they are no longer Grassmann variables. 

g) To summarize d, e and f , a variable living in the wrong countour is both a " Grassmann 
number" and a "variable", but it is not a "Grassmann variable". 

We claim that, particularly due to item c, the desired "rotation" properties can be 
accomodated both for the case of Grassmann numbers and Grassmann variables, although 
we need a lot less care in the former case than in the latter. 

When we are discussing symmetry transformations, we are discussing transformations of 
Grassmann numbers rather than Grassmann variables. Since the set of Grassmann numbers 
is a continuum, we can envision a "classical" fermionic field ip which varies continuouly and 
differentiably in space. However, we should not let the "continuity in x" confuse us. When 
we talk about "variable" we are claiming that, for any fixed x, the value of ip(x) can still 
"vary". This, of course, would be true in the context of path integral (and in this case 
we would be integrating ip(xi) over C±, ip(x 2 ) over C2, and so forth; or we could instead 
decide that they are all being integrated around the same contour C; but in either case, the 
statement about "contours" would turn them into "variables"). But for "classical" fermionic 
field, the value of ip(x = x ) is not assumed to "vary", which means it is not a variable. 
Since in case of Grassmann numbers we don't have to worry about "staying on a contour", 
we are free to rotate them. This, naturally, applies to "classical" fermionic field and allows 
us to do symmatry transformations in "classical" context. 

Now, of course, we would also like to be able to rotate Grassmann "variables" as well, 
and also we would like to think of an end-product as a " Grassmann variable" . For example, 
it is common place to replace the integral over 9\ and 82 with an integral over and 0, where 



and we would like to be able to think of the latter as something well defined in a literal 
sence, as opposed to mere convention. Let us start with "real valued" case, 



9 = 



V2 ' V2 



(31) 



91 + 92 
V2 




01 — 02 



(32) 
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In this case, we can simply agree that whenever different physical parameters are related 
by certain relations, the same applies to the contours they reside on. For example suppose 
we have a single spin 1/2 particle. In this case, we agree upon that the Grassann numbers 
related to spin-up and spin-down are "constrained" to the contours 6 = +z (a) and 6- z {a). 
Then the corresponding spins in +x and -x directions are constrained to contours given by 

9 +x (a) = ^=(ft+» + , 0-M = T^+» " M*)) (33) 

Now, due to the l/y/2 coefficient, the "new" contours obey the properties discussed in the 
previous section, which is why the "variables" associated with +x and —x are, in fact, 
Grassomann variables. The coefficients, in principle, could have been different. In this case, 
we would still have new contour, prescribed by the "wrong" coefficients, and we would still 
have new variables, but they would no longer be called "Grassmann variables". For example, 
if we assume that 

a) rj is a Grassmann variable 

b) X = 2V 

Then we would obtain 

J d X -X = ^ J dr 1 -r 1 = A (34) 

which would imply that x is n °t a Grassmann variable. It is also related to the issue of 
integrating of exponents. In particular, it is true that 

J drj ■ e 2r > = X - J d X ■ e x , where e A = 1 + A + + • • • (35) 

but, at the same time, 

d X ■ e x = 4 (36) 



/ 



which is why 



dr] ■ e 2r > = \ x 4 = 2 (37) 



Again, the number 4 is a consequence of the fact that \ is n °t a Grassmann variable. At 
the same time, \ ^ s still a "variable", and we are free to continue to use \ as l° n g as 
we remember the "new" integration properties. Furthermore, since x shares the same TV- 
dimensional space with r], it is still a Grassmann number. Thus, it continues to have all of 
the "algebraic" properties of dot and wedge products outside the integration; but the rules 
of integration are now different. 

Of course, what we have just said about i] and x is also true in the "conventional" 
presentation of Grassmann numbers. But in the "conventional" case the agreement as to 
when the integration properties hold and when they don't is merely formal. In our case, on 
the other hand, they are a consequence of the geometry of the specific contour. Thus, even if 
r] has never been mentioned, one can still know that the integral over x returns 4 rather than 
1 by simple inspection of the contour. This, in particular, means that the contour should 
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change the shape under linear transformations such as above. Since we can "stretch" straight 
line without changing its shape, this immediately tells us that the contour in question either 
has to have limitted length, or it has to have non-trivial curvature, or both. 



From what we have seen in the previous section, we already know that we need the 
contour to "bend around" for other purposes. In particular we wanted the properties of 
multiple integral to continue to hold true even if the same contour is selected to "constrain" 
all of the variables under discussion. In order for their wedge product not to return zero, 
we had to assume that the contour makes arbitrary many "turns". This, of course, has a 
side-benefit for the issue we are currently discussing. In particular, it means that we don't 
have to do "anything more" to assure that the "doubling" of a contour is distinct from the 
original contour. After all, the presence of different "turns" guarantees this to be the case. 

Finally, let us turn to complex case. As expected, this would require us to extend the 
space of Grassmann numbers from real valued to complex valued, which amounts to doubling 
the number of "real" dimensions: 



S = {nidi + ir 12 e u + r 2 ie 2 i + ir 22 e 21 H } (38) 

This creates a problem. As we have just mentioned, if we multiply a Grassmann variable by 
a "wrong" factor, it would no longer be a Grassmann variable. Now, % is an example of such 
a "wrong factor". For example, if x — iv then 

dx ■ X = - / dr}-r} (39) 



But, due to the difference in dimensionality, avoiding the factor of i would have imposed much 
more severe limitations on us than avoiding the factor of 2. Instead, we will deal with the 
issue by taking "complex conjugate" of the "differential" part, without taking corresponding 
conjugate on the "finite" part. Thus, if we assume that 

and if we also assume that 6>i and 9 2 are Grassmann variables, we obtain 

j S-e = i^l J de 1 ■ e 1 + jde 2 -e 2 = i (4i) 

On the other hand, 



/ 



d6 ■ 9 = / dd i ■ 0i + % -^r I d9 2 - 9 2 = (42) 



The reader should notice that the above two equations are just the opposite to what we 
would expect to get. The way we "fix" it is by introducing a third kind of product, o, 
defined as 

A o B = A ■ B (43) 
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Thus, we obtain 

J d9o9 = J d~9-9=l (44) 

J d9o9 = J d9-9 = (45) 

Since we are defining another kind of product for the purposes of " convention" , we might 
as well kill two birds with one stone and deal with one other "convention" that needs to be 
addressed. Going back to "real valued" case, it is usually assumed that 

d9 1 d9 2 9 1 9 2 = 1 (46) 

while in our case we have, by using quasi-anticommutativity of the dot-product between d9 2 
and 6*i in the "typical" scenario (that is, "typical" in a sense of special Manhattan metric 
discussed earlier), we obtain 

J d9 1 ■ d9 2 ■ (9 1 A 9 2 ) = - J d9 l ■ ^ • J d9 2 ■ 9 2 ^j = - J d9 1 ■9 1 = -l (47) 

We address the sign issue by claiming that 



A A B = B A A (48) 
Now, by substituting this into the definition of o, we obtain 

J (d9i ■ d9 2 ) o (9 1 A 9 2 ) = J d9[ ■ d9 2 ■ (#i A 9 2 ) = J d~9 2 ■ dB x ■ [9 X A 9 2 ) = 1 (49) 



In general, 



/ 



(d9 1 ■ . . . ■ 9 n ) o (9 1 A ■ • • A 9 n ) = J d9 n ■ . . . ■ d9 1 ■ {9 1 A • • • A 9 n ) = 1 (50) 



The left hand side of the above equation, indeed, agrees with commonly accepted convention. 
The price to pay for this, however, is the use of o which is not associative. 



5. Derivatives with respect to Grassmann variables 

Let us now turn to a much simpler issue and attempt to define the derivatives with respect to 
Grassmann coordinates. The only obstacle to overcome is the fact that we have to " divide" by 
"vectors". From convential point of view, whenever we compute df/d9, f is either constant 
in 9 or proportional to 9 (as always); thus, derivative is well defined in both cases. In the 
context of our model, however, we should be allowed to define some non-trivial function that 
can not be generated through wedge products; after all this is a simple consequence of the 
fact that we now have a set with well defined mathematical objects, rather than mere formal 
manipulation of symbols. Thus, it is possible for example, to define 

6>(A) = Ae 5 , /(Ae 5 ) = e 7 sinA, (51) 
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which means that the derivative will force us to evaluate e-jje^. The way to do it is to define 
the devision in terms of dot-product rather than wedge-product: 



^ = B ^9 ■ B = A (52) 

9 



in this case, we immediately obtain 



e 5 • (e 5 A e 7 ) = e 7 =>- — = e 5 A e 7 (53) 



In fact, the above definition of ratio is equivalent to 

A 9 ■ A 



9 9-9 



(54) 



where the right hand side is well defined due to 9 ■ 9 being a scalar. After all, by using the 
associativity of scalar multiplication on the first step, and associativity of the dot-product 
on the second step we obtain 

where on the last step we were using the assumption that 

k e C VA(k AA = AAk = k-A = A-k = kA) (56) 
Thus, if we have n Grassmann variables, 9 1 , - • • ,9 n living in the same contour 9(a), 

9 k = 9{a k ) (57) 

we can define 

F((Ti, • • • , a n ) — F((t(ti), • • • ,a(r n )) (58) 
we can formally define the partial derivative as 

OF 9\a k ) ■ (dF/da k ) 



d9 k 9\<j k ) ■ 9\a k ) 



(59) 



On the other hand, if the Grassmann variables live in the different contours (namely 9 k lives 
on its own contour 9 k (a) instead of the common contour 9(a)), one can simply rewrite the 
above expression while adding extra index k to all #-s: 

dF_ = 9\(a k ) ■ (dF/da k ) 
90 k 9\(a k ) ■ 9\(a k ) 

where one should remember that k is fixed, and therefore the "repetition" of index k does 
not imply summation. 

The above definitions would occasionally imply that the derivative has more components 
than we would expect. For example, we would normally expect the derivative of one vector 
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with respect to the other to return a scalar; yet, in some situations, including the example 
of Equation [5TJ we obtain a tensor instead: 



— = e 5 A e 7 cos A 



(61) 



If, on the other hand, we were to define g as 



0(A) = Ae 5 , g(Xe 5 ) = e t sinA, 



(62) 



then we would, in fact, obtain a scalar: 



cos A 



(63) 



Since in the actual physical situation we typically consider the functions that can be gen- 
erated through wedge products, we would not encounter the "surprising" scenario. But, in 
principle, it is possible to explore some non-algebraic modifications of the theory in which 
case we will begin to encounter them. 

Another thing that is important to mention is that, contrary to the situation with the 
integral, in case of derivative we do not need a contour in order for it to be well defined 
(although of course we were using the contour in the version of the derivative we described 
above). We can, instead, define the derivative in terms of a "gradient" with respect to entire 
space: 



where we have omitted the devision by • since we know that the latter is equal to 1. It 
should be noted that the "gradient" can be applied both to scalar, to vector and to higher 
rank tensor. If it is applied to scalar, we would get a vector as one might expect. Applying 
to vector will return a linear combination of a scalar and rank 2 tensor. Applying gradient 
to rank n tensor will return a linear combination of tensors of rank n — 1 and n + 1. 

In case of "taking the derivative" of path integral with respect to "source terms" in 
order to compute a correlation function, we can think of it as either "derivative with respect 
to a path" or "gradient with respect to entire space"; in either case we would get the same 
answer. From physics point of view, however, we understand that there is no such thing 
as "source term". Instead, we are dealing with one of the fermionic fields which, for the 
purposes of our own human convenience, we are separating from everything else. Now, since 
we postulate that physical fields are constrained to a contour, the same is true for the sources 
and sinks. Thus, from this perspective, we are taking the contour-based derivatives rather 
than gradients. But, again, this does not affect the actual answer we will arrive at. 

6. Harmonic oscillator in two Grassmann coordinates 

So far we have found a way to model grassmann numbers as mathematical objects. The next 
logical question is to do the physics in the Grassmann space we have just defined. Doing the 




(64) 



k 



15 



lattice quantum field theory on fermions is beyond the scope of this paper. Instead, we will 
look at the "first quantization" of the single particle living in two Grassmann coordinates. 
Since we have modeled Grassmann numbers as elements of a continuum, one would expect 
that the first quantization of physical system will look like a wave function, 

i/> = ip(e h e 2 ',t) (65) 

We will derive its evolution by means of path integral. Just like in bosonic case, in order 
to define path integral properly, one would need to discretize time. However, due to the 
anticommuting nature of Grassmann variables, the effects of the discretization will end up 
being much more significant than they were in commuting case. After all, 9i(t Q ) "commutes" 
with itself and, at the same time, it "anticommutes" with 9 (to + St). In fact, the product 
9 (to) A9(to + St) is "large" despite the fact that 9i(t ) A9i(t ) = 0. This opens some room for 
discontinuities and might possibly imply that our result will depend on the discretization of 
our choice. The study of different kinds of discretizations is beyond the scope of the paper. 
For our purposes, we will simply split time into even intervals of length St. Being inspired 
by bosonic path integral, we will postulate the "inductive" evolution of ip according to 

^(9i, 9 2 ; t = t + St)= [ (d9[ ■ d9' 2 ) o ( e -^MAA) A ^ ^. t = to) ) ( 66 ) 



where the exponentiation is defined in terms of wedge-products: 

A A A A A A A A A A 

e A = l + A+ — — + + ■ ■ ■ (67) 

2 6 

Let us now start at t = and see how it evolves. By substituting t = into Equation 
we obtain 



^j(9 u 9 2 ; t = St) = J (d9[ ■ d9' 2 ) o ( e -im,0'2,6x,9 2 ) A ^ 0/ . t = 0)) (68) 
If we now substitute t = St into Equation [661 we obtain 

V>(0i, 9 2 ; t = 2St) = J (d9[ ■ d9' 2 ) o (e-im&M*) A ^(e[, 9' 2 ; t = St)) (69) 
By substituting Equation [68] into Equation [69], we obtain 

V>(0i, 9 2 ; t = 2St) = [ (d9[ ■ d9' 2 ) o ( e -im^M) A 



A I / (d9'{ ■ d9'i) o (e-tmMM) A ^ t = o)) ) ) (70) 



Now, as we recall from our previous discussion, due to the fact that we have three different 
products, we can not simply pull e* £ ^ 1 ' e ' 2 ' 6,1 ' e ' 2 ^ inside the integral over 9'[ o 9' 2 and expect 
exactly the same result. At the same time, we have also previously argued that, in light of 
the special Manhattan metric, we will get approximately the same result, where the level of 
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approximation gets better and better as dimensionality of Grassmann space gets higher and 
higher. Thus, assuming that approximation, we rewrite Equation [7D] as 

A; * = 2£t) = J {do'l ■ de' 2 ' ■ d0 'i ■ do' 2 )o (71) 

The other suddle point is combining the exponentials. In general non-commuting case, the 
product of exponentials is given by Baker-Campbell-Hausdorff formula 

e A e B = e A + B + \[A,B]+- ( 72 ) 

However, in light of the fact that a typical Lagrangian contains only the products of even 
number of Grassmann variables, we can assume that it commutes with everything, 

C A A = A A C (73) 

If we make that assumption, Equation [72] implies that the product of exponents is just the 
usual exponent of the sum! Thus, 

By substituting Equation [7H into Equation [7H we obtain 

V(0iA;t = 2«) = y {dei-de^de[-de' 2 )o(e~ iC ^> e '^^ o)) (75) 

By simple induction argument, it is easy to see that, for arbitrary k, we will have 

V>(0i, t = = y ■ ^ 2 fe ■ • • • • d^u • d9 21 )o (76) 

( e -i£(flii,feiA.ft.)--i£(tfi*.fl afc ,ei,*-i.e2,*-i) a V(0ifcA*;* = 0)) 

This, indeed, shows that we obtain the path integral, as expected. 

Let us now attempt to tackle the problem of harmonic oscillator in fermionic coordinates. 
The first step is to write down a "classical" action. Concider a particle living in two- 
dimensional Grassmann space, and moving along a trajectory (#i(t)A(^))- I n a continuum 
setting, we would like the action to be 

S = J dtiwJ9x{t) A6 2 (t) -muO^t) A0 2 (t)) (77) 

Now, the definition of path integral we have provided is based on continuum space and 
discrete time. Thus, while leaving 9i and 9 2 continuous, we will discretize t. Thus, we will 
assume that 

t = k5t, ke {(),■■• ,n} (78) 
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We will replace "time derivative" at t according to 

^ e 2 ((k + m) -e 2 (kM) 

St 

In light of Grassmannian nature of parameters, it is conceivable that we would get different 
results depending on whether we chose "right derivative" or "left derivative". In order to 
be as general as possible, we will "split" constant terms between slice t and t + St. Thus, 
if we put "most" of the constant terms to t we would obtain "almost-right" derivative, 
while puting "most" of the constant terms to t + St would give us "almost-left" derivative; 
splitting constant terms "evenly" would give us "middle derivative". Thus, we will replace 
"same-time" terms according to 

77101 (*) ->■ 771101 (Arft) + m 2 0i((Jfe + !)<**) (80) 

mw6i(t) A 2 (f) ->■ mzuOxinSt) A 6 2 {nSt) (81) 
Where we have replaced m with three parameters m 1; m 2 and m 3 , given by 

m 3 = mi + m 2 (82) 

Intuitively, m 3 is the replacement of m, while rrt\ and m 2 are the two "parts" of m we 
"split" it as. Thus, if we want to have "left derivative" we would set = m 2 < mi = m 3 , 
for right derivative we would set = m\ < m 2 = m 3 , and for middle derivative we set 
< 777 1 = m 2 = 777 3 /2. By making these substitutions, discretized action (for general mi, 
m 2 and m 3 ) becomes 

71—1 

5 = ^(mi^i(Mi) A e 2 ((k + l)<5t) - m^kSt) A 6 2 (kSt)+ (83) 

fc=0 

+ m 2 0i((£; + l)<St) A e 2 ((k + l)<Ji) - m 2 0i((£; + l)St) A 6> 2 (fc5t) - m 3 w5t O^kSt) A 2 (ifc<ft)) 

We notice that the of the form 6\{{k + l)St) A 9 2 ((k + l)<5t) for /c < n — 1 match the terms 
of the form $i(k) A 2 (k) for A; — >■ k + 1. By combining such terms for 1 < /c < n — 1 while 
taking separate care regarding the "surface terms" (k — and k = n), we obtain 

S = m^nSt) A 6 2 (nSt) - m^i(0) A 2 (O) - m 3 w5t ^(0) A 2 (O)- 

rt-l 

- m^e^nSt) A 2 (n<J*) + ^(m^^Mt) A 2 ((Jfe + l)St)- (84) 

fc=i 

- 771201 ((A; + 1)5*) A 2 (Ht) + (m 2 - mi - m 3 uSt)9 1 (kSt) A 9 2 {kSt)) 
The above expression can be rewritten as 



fc=0 
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where we define Ck according to 

C = -(mi + m 3 uSt)9i(0) A 2 (O) 

C k = (mi6i(k8t) A 9 2 ((k + l)5t) - m 2 6 l ((k + l)5t) A 9 2 (kSt)+ 

+ (m 2 - mi - m 3 u8t)9i(kSt) A 9 2 {k8t)) , 1 < k < n - 1 (85) 

£„, = (m 2 — m 3 o;)0i(?7,<5t) A 9 2 (nSt) 

Let us now rewrite it in a notation of Equation [66j Thus, we will denote 6i(kSt) and 9 2 (kSt) 
by 0^ and 2 , while denoting 8i((k + l)St) and 9 2 ((k + l)5t) by Q\ and 2 . Thus, we rewrite 
C as 

£(0i, 2 , 6 U 6 2 ) = -(mi + m 3 u;<ft)#i A 6 2 , t = 

£(0£, 2 , ls 2 ) = (mi^i A 2 - m 2 0i A 9' 2 + (86) 

+ (m 2 - mi - m 3 u5t)6[ A8' 2 ) , 5t < t < (n - l)6t 

C(6[, 2 , 0i, 2 ) = (ma - m 3 ou)6[ A 2 , t = t n 

Let us now take the exponent of the above. As before, we will define the exponent in terms 
of a series of wedge products. This means that the expression terminates after linear terms. 
From straightforward calculation, we obtain 

e -ic = 1 _ + mqufyOi A 2 , t = 

- i(m 2 - mi- m 3 u)9[ A 9' 2 + mim 2 9[ A 2 A X A 9' 2 
e -ic n _ i _ _ m^s^nSt) A 9 2 (nSt) 

Let us now compute the evolution of ip. Suppose ip was arbitrary at t — (that is, it can 
be somthing completely non-linear). Now, any 0-dependence of ip at t = will become 0'- 
dependence as far as t = 5t is concerned. That dependence will disappear via the integration 
over 0'. On the other hand, new dependence will be created; namely, the dependence on 
un-primed 0-s. That dependence, however, will be strictly linear, since 0-s will enter linearly 
into Equation [HTl Thus, even if ip(9) "starts out" non-linear at t = 0, it will "collapse" into 
linear function by t = St. Therefore, if we are concerned about t > St, we can write down 
■0(01, 9 2 ) in a form of a general linear function: 

0(0i, 2 ; t) = 9Ji(t) + 9 2 f 2 (t) + 0i A 9 2 f 3 (t) + U(t) (88) 

The time evolution of ip, therefore, is completely described through the time evolution of fi, 
f 2 , f 3 and / 4 . Now, if we substitute Equations 1571 and into Equation I6"61 we obtain that, 
for St < t < (n — l)St, the dynamics of ijj is given by 

^(d u 2 ; t) = J (d9[ ■ d9' 2 ) o ((1 - imi9[ A 2 + zm 2 0i A 9' 2 - 

- i(m 2 - mi- m 3 u)9' 1 A 9' 2 + m x m 2 9' x A 2 A 9 X A 2 )A (89) 
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A (9'Mt - St) + 9' 2 f 2 (t - St) + 0[ A d'J 3 (t - St) + U{t - St))) 
By evaluating the above integral, we obtain 

i/>{6i, 9 2 ; t) = -im 2 6 l }\{t - St) + %m x Q 2 j 2 {t - St)- 

— m\m 2 Q\ A 62f4.it — St) + fs(t — St) + i(m 3 uj + mi — m 2 )f±(t — St) 
By comparing Equation [88] to Equation [90], we obtain 



(90) 



/ fi(t) \ 

hit) 
hit) 
\ hit) J 



( -im 2 

inn 

—m\vri2 

\ 1 i(mi — rri2 + m^uiSt) J 



\ ( hit -st) \ 
hit -st) 
hit - st) 
V hit - st) ) 



(91) 



Thus, naively we would expect that e lESt coincides with the eigenvalues of the above tensor. 
Now, the eigenvalues associated with the top left block are 

(92) 



c\ = -yjmiin 2 , c 2 = A/mim 2 
while the eigenvalues associated with bottom right block are 

i{nii — m 2 + m 3 uSt) — \/—(mi — m 2 + m^uSt) 2 — 4mim 2 



c 3 



C4 



i(mi — m 2 + m^uSt) + y/—{rrii — m 2 + m 3 uiSt) 2 — 4mim 2 



(93) 



(94) 



Now, we assume that m-s as well as u are real, it is easy to see that the expression under 
square root for C3 and C4 is negative, which makes C3 and C4 purely imaginary. At the same 
time, ci and c 2 are purely real. Thus, if we stick to the view that 



we would be forced to conclude that 

-Ei, 2 



7T 



7T 



3,4 



2w5t 



(95) 



(96) 



This of course has several major problems. Just to name a few, we have four eigenvalues 
instead of two, in contradiction to Fermi exclusion principle. In fact, two of the eigenstates 
would have twice larger energy level than the other two, while in a simple fermion system we 
would expect the two energies ("vacuum" and "one particle") to have the same absolute value 
and opposite sign. Finally, the energies are simply too high. The fact that the "periods" of 
the "oscillation" processes match either St or 2St would not account for any "more obviously" 
sinusoidal processes seen in the lab. We avoid the above issues by proposing more clever 
way of defining energies: 

a) Assume that u is imaginary. In fact, set u = iE (for now, E will simply be a letter 
in the alphabet but we will soon show this in fact represents energy). 
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b) Factor out the overall complex phase in the eigenvalues and look at the deviations in 
phase due to u = iE. 

c) Try to see whether or not the overall factors we have "ignored" in part b might lead 
to some problems, such as interference between states and so forth. 

This immediately tells us that the states associated with eigenvalues c\ and c 2 have zero 
energy, since u has no impact on their values: 

E x = ? E 2 = ? (97) 

At first, this might appear troubling since we know that fermionic vacuum has negative 
energy rather than 0. Upon further look, however, one realizes that due to the block- 
diagonal form of the matrix, these two states are "isolated" from the states associated with 
C3 and C4. Furthermore, since cu-term originally came from coupling, the fact that states 1 
and 2 lack that term shows that they are neither coupled to themselves nor to anything else. 
Thus, even if we were to perturb our system in some way, that perturbation will amount to 
the "additional interaction" between states 3 and 4, while states 1 and 2 will continue to 
stay "isolated". This simultaneously expins why we "observe" only two states and why we 
"don't observe" zero-energy that is apparently attributable to states 1 and 2. 

As far as states 3 and 4 are concerned, the evaluation of Equations [93] and [94] up to 0(r) 
tells us that 

im 3 u5t i(m\ — m 2 )m 3 uj5t , . 

c 3 = —im,2 H ; : 98 

2 2(m 1 +m 2 ) V ; 

im 3 u5t i(mi - m 2 )m 3 uj5t 
c 4 = imi + — _ — + — _ . 99 

2 2(m 1 +m 2 ) 
Without loss of generality, we can substitute 

m 3 = nil + m 2 (100) 

After all, we still assume that m-s are real, while any mismatch in amplitude of 777.3 can be 
absorbed into the amplitude of u. Thus, we go ahead and substitute the Equation 11001 into 
the Equations EH and [HH1 to obtain 

C3 = — 2777,2(1 — oj5t) , C4 = im\{l + uiSt) (101) 

which, upon substitution uj = iE becomes 

c 3 = -im 2 (l - iESt) , c 4 = zmi(l + iESt) (102) 

This looks a lot more promising. If we were to ignore the overall coefficients and just look 
at the terms inside the bracket, we find that, indeed, we have a "very slow" time evolution 
that would produce sinusoidal process. Furthermore, we have E comming with the same 
magnitude and opposite signs, just as expected from vacuum and one-particle states. Now, 
if the overall factors outside of the bracket were the same for C3 and C4, we would have been 
able to just normalize them away. Due to the fact that they are not the same, however, poses 
a possibility of interference between c 3 and c 4 that would "disclose" the unwanted "rapid 
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evolution". However, in case of m\ = m 2 the interference will simply amout to relative 
sign between C3 and C4 alternating between 1 and —1 in aforegiven order; its effect can be 
reasonably expected to average out to something trivial. If we assume that m% 7^ m 2 things 
might get a little more interesting since the state with larger m will dominate over the state 
with smaller m; in fact, the effect of the latter will quickly approach zero in comparison to 
the former. We will, therefore, have to assume that 

m i — m 2 (103) 

in order to avoid the above situation. This setting is most natural anyway since it respects 
time reversal symmetry. If we now substitute Equation 11031 along with 

m = 1713 , u) — iE (104) 

into the Equations [98] and [99j we obtain 

—Em5t — \/m 2 E 2 (St) 2 — m 2 
c 3 = (105) 



-EmSt + Jm 2 E 2 (St) 2 - m 2 . . 

c 4 = (106) 

Now, since St is "very small", the expression under square root is negative, which means 
that square root itself is imaginary. Now, since the square root contains m, which is finite, 
while the term outside of square root has St in it, the square root dominates. Therefore, due 
to imaginary nature of square root, it makes sense to factor out im rather than m. This 
implies 

c 3 = -™(y/l-E*(6t)*-iE6t) (107) 
C4 = ^ (^l - E 2 (Stf + iESt) (108) 

By re-expressing this as 



= Wl-W / _ iESt 

2 V y/i - E 2 (sty 1 



imJl - E 2 (St) 2 f iESt \ , . 

c 4 = — ^ 1 + ; 110 

2 V ^l-E 2 (St) 2 J V ^ 

which shows us that the higher order corrections due to square root do not imply the dif- 
ference between the absolute values of positive and negative energy states. However, apart 
from this curiocity, we areonly computing the subsequent exponential evolution up to finite 
order. After all, due to the higher order deviations between 1 — iESt and e~ lESt , doing higher 
order calculations would be increasingly more complicated; in this respect the situation is 
no different from bosons in discrete time. Up to finite terms, we see that the two eigenstates 
evolve as 

/ • \ t/ St 
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■ \ t/St 



(112) 



It is easy to see that the corresponding eigenstates can be expressed as 



v 3 (0) = C 



( o 



—i 

\ 



v 4 (0) = D 



J 



/0 


i 

\ 



(113) 
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The above is not to be confused with spinor. As was stated earlier, the only "physical" 
components are the third and fourth component of the above column, whereas the first two 
components are spurious. What we are describing corresponds to one spin-up particle state 
and vacuum state. The "spinor" would include four copies of the above (corresponding to 
spin- up particle, spin-down particle, spin- up antiparticle and spin-down antiparticle); thus, 
it would be having 16 components, 8 of which physical and the other 8 spurious. In this 
paper, for the sake of simplicity, we restrict ourselves only to one spinor component, which 
is why we get only 4 states (2 physical, 2 spurious) instead. Any similarity between what 
we have and spinor notation is purely coincidental. 

Now, if we substitute Equations 11111 and 11121 into Equation 11131 we obtain 
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This means that corresponding wave functions evolve as 
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For the completeness purposes, let us now compute "spurious states" ipi and ip 2 - It is easy 
to see that their eigenvalues are C\ = —m and c 2 = m, which correspond to 



vi(0) = A 
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This means that these states evolve as 
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The lack of any exponential term in the coefficient is closely related to the fact that these 
states are spurious. Finally, they correspond to the wave functions ipi an d 4>2 evolving as 

M&i, #2; t) = Ai-mf^iiet + 2 ) , ^i, 2 ; t) = Bm t ' St {9 l + i9 2 ) (119) 

The general wave function, therefore, evolves according to 

V>(0i, 2 ; *) = Ai-mY^ite, + 2 ) + Bm t ' 5t {9 l + ^ 2 )+ (120) 

The terms with coefficients A and B are spurious, the term with coefficient C represents 
one-particle state and the term with coefficient D represents vacuum state. The spurious 
states have zero energies, the one-particle state has energy E and the vacuum state has 
energy —E. 



7. Coupling between bosonic and fermionic coordinates 

Let us now discuss the situation when we have three coordinates: two Grassmannian and 
one real, and see what happens. The wave function will evolve according to 



o e 



if>(6 1 ,6 2 ,x,t) = J \(dx'd9' 1 -d9 2 )o (121) 
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Now, let {ip n (x)} be an eigenfunction of Schrodinger's operator; that is, 

1 d 2 1p n 



2m dx 2 



+ V^ n = E n ^ n (122) 



The above step is pure mathematics: since we have extra anticommuting coordinates, we can 
not claim that Schrodinger's equation holds, unless we have demonstrated it through path 
integral. But, from the point of view of differential equations, the above set of eigenfunctions 
continues to exist. Since that set of eigenfunctions "happen" to ber complete, the general 
analytical function ip{9i, 2 , x, t) can be written as 

^(9 1 ,9 2 ,X,t) = Y, d ^n{x)f ln (t) + Y^ d ^n{x)hn{t) +^«lA 9 2 f 3n (t) + Y,Mx)Un(t) 

n n n n 

(123) 

Let us now evaluate Equaton ll21l by substituting Equation ll23l for if}(0i, 62, x, t) and Equation 
[ST] for e - lC ( e 'i< e 2' 9 ^M , By evaluating the integral over 9[ ■ 9' 2 we obtain 



ip(9i,e 2 ,x,t) 



£ / dx'e~ i < ! ^- V ^(M^fsn(t-5t)+ 
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+ %m x Q 2 ^ n {x')f 2n {t - St) - im 2 9 1 ^ n (x')f ln (t - St)- (124) 
- i(m 2 - mi - m 3 u)ip n (x')f An (t - St) - m 1 m 2 6 1 A 8 2 ip n f An (t - St) + ill (x')i/; n (x')f An (t - St)) 
Now, in light of the fact that {ipk] is complete, we can express U(x)ip n (x) as 



U(x)tp n (x) = ^C kn i) k {x) 



(125) 



where 



Ckn= / U(x)ip*(x)ijj(x)dx 



(126) 



Furthermore, from our knowiedge of correlation between bosonic path integral and Schrodinger's 
equation, we know that 



/ 



dx e 



st 



m£} (x — x ) 
2(5t)' 2 



-V(x>)) 



x 



2nSt 



e- lE " St 6Jx) 



im 



By applying Equation 11251 and 11271 to Equation 11241 we obtain 
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+ im x Q 2 ^n{x)hn{t - St) - im 2 dx^n{x)f Xn {t - St)+ 
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Let us first consider the case where bosons and fermions are decoupled; that is, we will assume 
that Ckn = 0. Now, by comparing the Equation 11281 to Equation 11231 we see that the step- 
by-step evolution of terms with "fixed" index n in this section is similar to the step- by step 
process of "everything" in the previous section (this similarity is "broken" via C nk terms, 
which is why we are first considering the case of these terms being zero). Therefore we can 
"copy" the Equation 11201 and insert the indeces n as well as a summation over them. Apart 
from that, we also have to be careful with coefficients that we now hae which we did not 
have before. As far as the ftnSt/im) 1 / 2 coefficient goes, that "extra multiplication" occurs 
every time we pass the interval of St. Thus, we have to take that coefficient to the power 
of t/St, which would produce (27iSt/imy^ 25t . As far as e~ lEnSt , in the expontial we would 
produce a sum of "very large" (t/St) number of copies of "very small" numbers (—iE n St), 
which would produce a "finite" number —iE n t in the exponential; that is, we get an extra 
coefficient of e~ lEnt . Thus, the "non-interacting" fermionic field (with Ckn = 0) will produce 
a copy of Equation 11201 with the insertions of summation over n as well as the coefficients 
just discussed: 
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By combining the e~' lEnt on the outside with the e~ tEt and e lEt on the inside, we obtain 

-^-J ^ U(-m) t / St (i0 1 + 9 2 ) + Bm t / 5t {9 1 + i9 2 )+ (130) 

This tells us that we obtain eigenstates with energies E n + E and E n — E 1 . The energy 
.E^ > comes from bosonic part, whereas subtraction or addition of E comes from negative 
energy "fermionic vacuum" and positive energy "fermionic one particle state" respectively 
(the bosonic vacuum energy, on the other hand, is absorbed into E n ). 



8. Multidimensional harmonic oscillator and multiparticle states 

So far we have worked out the "harmonic oscillator" in two Grassmann coordinates (which 
is the lowest number of coordinates where we could include "mass term" into). Now,the 
simple harmonic oscillator generalizes to multidimensional one which serves as a model of 
second quantization of bosonic fields. Multidimensional bosonic system continues to have 
wave function representation, except that it is now a wave fuction in N coordinates, where 
N is the number of lattice points. Suppose, for example, we have one-dimensional chain of 
points. Then their potential is given by 

l ' = 2?s3*D*»'-« ,+I r£tf (131) 

^ > k=l k=l 

which, upon simple agebra becomes 
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By diagonalizing the coefficients, we see that in "rotated" coordinates V is given by 

k=\ 

This means that the wave function is given by 

D 



V = £ ™- , where <ff k = £ (133) 
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where ip u - n (x) is the n-th eigenfunction of simple harmonic oscillator of frequency u. If we 
like, we can rotate the wave function back into un-prime frame, thus obtaining 

Ad) = ( C ^W ,n D Y[^ k ;nJ J^B^A J (135) 

ni,- ,n/j ^ fc=l ^ i ' ' 

where 

A ik B kj = 5) (136) 

Now, we "independently know" that the eigenstates of non-interacting field are fixed mo- 
mentum states. Thus, we can guess that (f>' k represents a particle moving with a "fixed 
momentum" p k , whereas <p k represents the particle of fixed position x k . This means that 

A kl » e ipkXl , B kl « e"** 3 " (137) 

up to some trivial power of 2ir. This essentially tells us that, even if we are "believers" in 
position, we can still "take seriously" wave function in a momentum-based Fock space. After 
all, Fourier transformation is merely a "rotation" inside the multidimensional space where 
our wave function lives. Since we agree that in "few" dimensions we are free to "rotate" 
wave function, the same is true for many dimensions. Thus, a general bosonic state 

w = £ E E ^ 1 ,-™,-nj<) ni ---(aL) n i°) ( 138 ) 
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We can now generalize it to fermionic case. In Section 6 we have established that the 
fermionic harmonic oscillator can be viewed as a wave function in two Grassmann coordi- 
nates. This, in conjunction with what we have just said in the present section, suggests 
that multiparticle fermionic system can be viewed as multidimensional Harmonic oscillator 
which can likewise be modeled in terms of wave function. Furthermore, in Section 4 we 
have established that we are "allowed" to rotate Grassmann variables and, provided that 
the roation is unitary, we can continue to think of end product as Grassmann variables, as 
well. Therefore, we can interpret Fourier transformation as that very rotation (after all, 
we have already agree that we can think of it as "rotation" in bosonic case and all of the 
obstacles that prevent us from generalizing it to fermions have been removed in prevous 
sections). Now, a general fermionic wave function in two Grassmann coordiantes has two 
physical states takes the form 

^(0i,0 2 ;t) = + / 2 (*)0 2 + / 3 (*)0i a e 2 + f 4 (t) (140) 
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where fi(t) and f2{t) are spurious parameters, while fs(t) and j\(t) are physical. Now, 
(61,62) represents a single spinor component in a single-point lattice. Upon "rotation" 
into momentum space, we can say that (61,62) represents a spin-up particle with fixed 
momentum. In order to account for spin-down particle as well as antiparticles of both spins 
of that given momentum, we need 8 coordinates, (6\,--- ,6s)- Finally, since we have D 
values of momentum; after all, we have D vaues of position comming from the lattice, so 
their "rotation" in a form of momentum should also have D values. Thus, the total number 
of 6-s is 8D. Thus, the general wave function is represented by 

SD 

^(6 l7 --- ,6 8D ;t) = J2 E C kl ... k J kl A---A6 kn (141) 

n=0 l<fci<-<fc„<8£) 

Now, if {2j + 1, 2j + 2} C {ki, ■ ■ ■ , k n }, this means that the "slot" number j is occopied 
by something "isomorphic" to 6\ A 62 which happens to be a linear combination of two 
"physical" states (a physical particle and physical vacuum). Likewise, if {2j + l,2j + 2} fl 
{ki, ■ ■ ■ , k n } = 0, then the j-th "slot" is occupied by a state isomorphic to 1, which another 
linear combination of the above two physical states, and therefore also physical. On the 
other hand, if 2j + 1 is an element of {k±, ■ ■ ■ , k n } but not 2j + 2, the j-th slot will be filled 
by something isomorphic to 6±, which is spurious. If the opposite of these two statmements 
is true then the j-th slot will be filled by something isomorphic to 6 2 which is, likewise, 
spurious. 

We can physically interpret the above situation in the following way. As usual, the 
vacuum energy is —1/2, while the energy of each particle is +1. But, on top of this, we also 
have "spurious particles" whose energy is +1/2. Thus, if the vacuum is occupied by the 
ordinary particle, the energy will be —1/2 + 1 = +1/2. If, on the other hand, it is occupied 
by spurious particle, the energy will be —1/2 + 1/2 = 0. Spurious particle is also a subject 
to Fermi exclusion principle. Thus, it prevents any other particle from filling a given slot; 
but at the same time, it does not interact with anything else either. This has several effects. 
First of all, it is possible that a given latticce point is occupied by spurious particle, but not 
spurious antiparticle. As a result, an antiparticle will be able to pass by that point, put a 
particle won't. Or, if a given lattice point is occupied by both spurious particle and spurious 
antiparticle, of both spins, then a fermion can not pass that lattice point. However, that 
lattice point can still be occupied by a boson. 

In addition to this, there is no interaction between different states whose spurious par- 
ticles fill different spots. As a result, these states might be macroscopically distinct from 
each other, which raises a question why don't they interact gravitationally. This question, 
however, can be asked in the context of any other quantum mechanical model involving de- 
coherence. The ultimate answer to such question is some form of collapse mechanism. It is 
quite possible that we might take advantage of the fact that Grassmann numbers are defined 
outside path integral i order to go beyond the realm of analytic functions in proposing such 
mechanism. 
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9. Conclusion and outlook 

In this paper we have defined Grassmann numbers in such a way that they are well defined 
outside path integral, which in principle allows for the realistic interpretations of fermionic 
fields. It turns out that if we properly define dot and wedge product, as well as "special 
manhatten" metric, then we would be able to expect the "limit of the sum" to approxi- 
mate Grassmann integral up to some factor, provided that the path we are integrating over 
has some very long series of mostly non-overlapping "turns" in very large number of non- 
commuting dimensions. Furthermore, we have worked out the behavior of a toy system in 
two fermionic dimensions. The discrete nature of time ends up playing a crucial role in 
the dynamics. We assume the discreteness scale to be St; this means that any given point 
of time is an integer factor of the latter. At slice t — 0, wave function can be arbitrary. 
Then, by the time t — St it "collapses" to something linear. After that, the linear function 
undergoes some periodic finite changes every time interval St. These changes, however, fit a 
trivial pattern and, therefore, deemed to be unphysical (although we would have been facing 
extreme discontinuity of the theory were we to take them seriously). However, this "quickly 
changing" function is being multiplied by "slowly changing" harmonic, and that harmonic 
represents the physical evolution of the state. It turns out, however, that only the constant 
term and 9 1 A 9 2 term are multiplied by corresponding harmonics, while the 9\ and 9 2 terms 
are not. Thus, the former two terms are interpretted as physical while the latter two terms 
are interpretted as spurious. It turns out that both vacuum state and one particle state are 
linear combination of the two "physical" components we just mentioned, with coefficients 
specified in this paper. 

One surprising result of our model is the predction of "spurious states". The root of 
the issue is the fact that we need two variables (6\ and 9 2 ) in order to obtain two states 
(vacuum state and one particle state). However, once we have these two variables we have 
four possible linear functions rather than two: 9 1: 9 2 , 9 X A 9 2 and 1. Thus, only two of these 
four degress of freedom are physical, while the other two are suprious. It turns out that 
the two physical degrees of freeom are 1 and 9\ A 9 2 , both of which correspond to linear 
combination of vacuum and one particle states with the coefficients given in the paper. The 
other two states are "spurious" and they don't interact with any of the physical states. 
As a result, it is possilbe for "spurious partcle" to permanently occupy one of the lattice 
points, thus not allowing the corresponding fermions to be created there, while the bosonic 
fields would continue to pass by that lattice point. This would potentially produce some 
non trivial deviations from lattice field theory that are worth investigating. On the other 
hand, the "exact match" with conventional theory is still possible: namely, in a scenario 
devoit of spurious particles. Since spurious particles do not interact, their absence "at the 
beginning" implies that they will never appear at all. Thus, if we are desperate to match 
the results of conventional physics we can simply postulate the absence of spurious particles 
as initial conditions. Such initial conditions, however, are quite unnatural and in my view 
it would be interesting to allow the spurious particles to continue to exist and explore the 
various possible outcomes that might be produced. In fact, it might even be interesting to 
ask whether or not doubling problem will persist in such a radically different discretization 
of fermions (currently I am of no opinion on this issue, so it is just a blind guess). 
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One of the originally intended benefits of our framework is that we gain an ability to 
view Grassmann variables in a literal physical sense, which might potntially provide some 
extra insight to interpretation of quantum mechanics. However, our theory inherits some of 
the limitations that we expect from ordinary fermionic theory. In "conventional" case, the 
limitation on the number of fermionic states corresponds to the lack of continuous spectrum 
of Grassmann variables. In our case, we do have continuous spectrum but, instead, our 
dynamics is limitted to producing only linear functions on that spectrum. In conventional 
case, we can't be speaking of "smoothly varying" fermionic field because we would require to 
have continuous spectrum in order to define the latter. In our case, we do have continuous 
spectrum but the multitude of "turns" on Grassmann space create their own problem. As 
soon as the two Grassmann variables on a given contour are separated by more than one 
"turn", their distance would be expected to be in the same range regardless of whther they 
are separated by only one turn or a million of turns; this is particularly due to the Manhattan 
metric that we have defined. Thus, in order to truly speak of Grassmann varriables being 
"very close" to each other, we have to assume that they are within the same "very small" 
part of the "turn" of the contour. In this case, however, most of the properties of Grassmann 
integration would disappear and it is questionable whether or not the physical implications 
of this would be of any interest. 

On a good side, however, it is still possible to use the extra degrees of freedom offered by 
the continuous spectrum to come up with some physical applications. While it is true that on 
a scale of several different turns of the contour we can not speak of "continuous functions" 
in a classical sense, we can still find other creative ways of introducing extra degrees of 
freedom. For example, we can consider a possibility of function deviating from linearity in 
some way that correlates either with the shape of the contour or with the deviation of some 
other function from linearity. Now, the specific dynamics that we have proposed in this 
paper (which was derived from the " naive" use of path integral) would kill any non-linearity 
during the very first interval St. It is, however, possible to explore some modifications of 
the definition of "time derivative", particularly the ones that can not be generated by the 
products, and see whether or not it is possible to retain some of the non-linearities. Of 
course, if these modifications are assumed to be small, then the non-linearities should be 
"almost killed" within "very short but not as short" period of time. This would, however, 
impy the possibility that during the big bang these non-linearities were still present in a 
non-negligible form. This would imply that Fermi exclusion principle did not work back at 
the time of big bang, which might lead to some novel cosmological theories. 

Appart from that, it is possible to attempt to apply this paper to the following problems 
(although I personally feel it might be highly questionable due to the issue with "turns" 
mentioned earlier): 

I. Causal Set Theory It has been recently attempted to define quantum field the- 
ory on a causal set (see [B]). In order to define Lagrangian density, we needed to consider 
Lorentzian neighborhood of any given point, which, in discrete case, has to be finite. Lorentz 
covariance implies that any finite Lorentzian neighborhood fills up the vicinity of light cone 

1 A causal set is a model of Lorentz covariant discretization of spacetime manifold which was originally 
proposed by Rafael Sorkin (see [3]) 
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and, therefore, is highly non-local. Bringing back locality would require "cutting off' the 
Lorentzian neighborhood and, therefore, violating relativity. In order to stay loyal to Lorentz 
covariance, we decided to come up with "Lorentz covariant" criteria of "selecting preferred 
frame". Namely, for any field configuration, the frame where field varies the least is "pre- 
ferred". This requires us to look at a "fixed" field configuration outside of integral. We then 
repeat this "analysis" for every single field configuration and then "sum up" our results at 
the end. In order to generalize the above from bosons to fermions, we have to be able to 
ask ourselves "in which frame does fermionic field varies the least". This question can be 
asked and answered only if the values of fermionic field (that is, Grassmann numbers) are 
well defined outside of path integral. 

II. Pilot Wave Models There were recent attempts by Struyve and Westmann (see 
[9]) to come up with pilot wave model for bosonic fields, where the fields, themselves, are 
viewed as "beables". However, the fact that fermionic fields are Grassmannian makes it 
very difficult to come up with similar model for fermions (although an attempt to do so was 
made by Valentini in his thesis). Instead, fermions are typically modeled to be particles (see, 
for example, [7]). Havind defined Grassmann numbers might now allow us to explore the 
possibility of field beables for fermions. 

III. Lattice theories As one would gather from the body of the paper, one of the 
underlying assumptions is that the dimensionality of Grassmann space is "much larger" 
than the number of variables we are integrating over. However, in a real physical situation, 
we are integrating over fermionic fields at every single point in space. Thus, in order to stay 
consistent with this paper, we have to insist that the space is discrete and, furthermore claim 
that the number of discrete points is much smaller than the dimensinality of Grassmann 
space, which is still assumed to be finite! For future research, it might be interesting to 
attempt to relax the last assumption. In this case, we might possibly predict some non- 
trivial differences from the "typical" lattice models of fermions. It might, for example, 
be interesting to explore whether or not the "new" version of lattice model will still have 
doubling problem. At the moment, the author of this paper has no opinion one way or the 
other. 
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